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摘 要 本文 利 用 CFD 数值 模拟 软件 对 当量 直径 为 200 hm 的 不 同和 孔隙 率 的 又 排 圆 形 微 针 肋 热 沉 的 流动 和 传 热 特 性 进 
行 了 数值 模拟 . 模拟 结果 表明 ， 随 着 雷诺 数 的 增加 ， 针 肋 尾 部 出 现 不 同形 态 的 涡 . 不 同 孔 隙 率 的 通道 中 ,由 于 后 排 针 肋 的 
影响 ， 针 肪 尾 迹 区 的 涡 的 形态 、 大 小 出 现 不 同 。 对 不 同和 孔 际 率 的 微 针 肋 热 沉 ， 热流 密度 的 变化 对 热 沉 的 努 塞 尔 数 Nu 的 
影响 不 同 , 孔隙 率 较 小 时 热流 密度 对 热 沉 的 努 塞 尔 数 Nu 影响 不 大 , 孔隙 率 较 大 时 热流 密度 升 高 , 热 沉 的 努 塞 尔 数 Nw 升 高 . 
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Abstract In this paper, by utilizing CFD, numerical simulations and analyses were used to research 
the flow and heat transfer characteristic in different porosity staggered circular shaped micro pin-fin 
heat sinks whose hydraulic diameter is 200 um. The results indicate that, with the increase of 
Under 
different porosity, the size and shape of vortexes are different due to the influence of pin-fins at 


Reynolds number, vortexes with different structures are found in the pin fins wake zone. 


downstream. The influence on Nusselt number of the heat sinks caused by the change of heat flux is 
different under different porosity. There is barely influence on Nusselt number under smaller porosity 
when heat flux changes, while under larger porosity, the Nusselt number increases with the increase 


of heat flux. 
Key words numerical simulation; micro pin-fin; heat transfer characteristic; low distribution anal- 
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随 着 微 电 子 技术 的 快速 发 展 ， 微 流动 系统 在 高 
热流 密度 微型 咒 件 及 设备 的 散热 中 得 到 越 来 越 广泛 
的 应 用 . 作为 微 流动 系统 的 重要 组 成 部 分 , 微 通道 热 
沉 因 具 有 较 好 的 强化 传 热 效 果 ,得 到 了 国内 外 学 者 
的 广泛 关注 和 研究 ， 并 发 展 出 强化 传 热 效 果 更 好 的 
微 针 肋 热 沉 . Kosar 等 0 和 首先 开展 了 微 针 肋 热 沉 
流动 换 热 方 面 的 工作 ， 对 微 针 肋 热 沉 的 摩擦 系数 的 
计算 进行 了 研究 ,提出 了 适用 于 微 针 肋 热 沉 的 摩擦 
系数 关联 式 . 随后 , 用 模拟 和 实验 的 方法 , 对 不 同形 
状 、 不 同 布置 方式 的 微 针 肋 热 沉 的 传 热 与 压 降 进 行 
了 研究 , 指出 通过 合理 地 设计 微 针 肋 热 沉 , 可 以 在 保 
持 优 越 传 热 性 能 的 基础 上 有 效 减 小 流动 阻力 ， 进 而 


收 稿 日 期 2017-02-16; 修订 日 期 ，2017-07-20 
基金 项 目 : 国家 自然 科学 基金 (No.51576005) 


作者 简介 : 杨 宇 撒 (1991- ), 男 , 硕士 研究 生 , 主要 从 事 微 扩 度 强化 传 热 研究 . 


降低 泵 功 消耗 , 周明 正 等 站 和 夏 国 栋 等 @ 分 别 对 流 
体 模 擦 水滴 形 微 针 肋 热 沉 和 获 形 、 长 萎 形 微 针 肋 热 
沉 进行 了 实验 研究 , 研究 表明 ,在 相同 布置 方式 下 ， 
相对 于 贺 形 和 萎 形 针 肪 ， 水 滴 形 和 长 萎 形 针 肋 在 设 
计 上 避免 了 尾 涡 澡 留 区 的 形成 具有 较 好 的 传 热 效 
有 果 。 由 于 在 微 电 子 设备 的 实际 应 用 中 ， 屁 件 的 发 热 
量 会 随 工 作 条 件 的 变化 而 变化 ， 对 热 沉 的 散热 效果 
产生 影响 , 因此 一 些 研究 人 员 也 开始 研究 热流 密度 
变化 对 热 沉 换 热 性 能 的 影响 . Guan 等“- 引 对 不 同 
热流 密度 下 圆 形 针 肋 尾 部 旋涡 的 变化 ， 以 及 热 沉 中 
边界 层 厚度 的 变化 进行 了 研究 ， 发现 热流 密度 升 高 
会 使 针 肋 尾 迹 区 拉 长 ,边界 层 减 薄 ， 壁面 效应 减弱 ， 
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使 热 沉 对 流 换 热 能 力 增 强 。 


当前 ， 对 微 针 肋 热 沉 的 研究 多 集中 在 传 热 和 压 
降 特 性 方面 ， 对 不 同 排 布 方式 下 多 圆柱 绕 流 的 对 比 
研究 还 较 少 ; 同时， 热流 密度 变化 对 不 同 排 布 方 式 
的 微 针 肋 热 沉 的 影响 也 不 明确 。 因 此 ， 本 文 用 数值 
模拟 的 方法 , 在 不 同 热流 密度 下 , 对 当量 直径 为 200 
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1.2 控制 方程 

为 简化 计算 ,本 文 对 模型 做 出 如 下 假设 : 三 维 、 
稳 态 、 无 内 热源 层 流 流动 和 传 热 过 程 ， 工 质 颖 度 随 
温度 变化 按 分 段 线 性 变化 ; 忽略 辐射 和 热 损失 作用 . 
由 假设 可 知 ， 流 体 的 连续 性 方程 、 动 量 方程 和 能 基 
方程 可 简化 为 ; 


um, 孔隙 率 为 ==0.743 和 s=0.836 的 两 种 微 针 肋 热 Ce =0 (1) 
yi ~ 执 特 性 济 行 宗 
沉 中 针 肋 扰 流 流 场 和 传 热 特 性 进行 了 研究 ， Blpruawy) _ 9P) 9 E | 
O(peus 0 0 
1.1 物理 模型 ee (天 
考虑 到 两 种 微 针 肋 热 沉 模型 为 对 称 结构 ， 因 此 (3) 


选取 模型 的 一 半 作为 模拟 单元 , 模拟 基本 单元 如 图 1 各 ) + ( 星 : 
所 示 , 针 肋 正面 结构 如 图 2 所 示 . 针 肋 排 布 方式 如 表 Os Oa 
1 所 示 . 根据 不 同 的 针 肋 排 布 方式 , 两 个 圆 形 微 针 肋 。 固 休 区 域 能 量 方程 为 : 
热 沉 的 孔隙 率 分 别 为 。 = 0.743 和 = = 0.836, 将 两 个 
微 针 肪 热 沉 分 别 命名 为 SCP-0.743 和 SCP-0.836( 其 Oz; ( 
中 S 代表 又 排 布置 ，C 代表 圆 形 , P 代表 针 肋 , 数字 
代表 孔隙 率 ). 固体 基板 高 万 =0.15 mm， 加 热 面 长 
一 工 =10 mm, 宽 WW=1.76 mm, 针 肋 高 82=200 um, 前 导热 系数 ,pr 流体 运动 和 灶 度 ，》s 为 固体 导热 系数 ， 
一 ”后 各 有 3 mm 的 稳 流 区 , 不 进行 加 热 。 1.3 边界 条 件 

ba 假设 流体 与 固体 接触 面 没 有 速度 滑 移 . 边界 条 


Ou; 2 
Oj 


oT 
As 区 ) =0 (4) 


其 中 , pr 为 流体 密度 ,cp,t 为 流体 比热容 ,Xi 为 流体 


F- 件 设置 为 
入 DI: 
T=0, uu=un, v=0, w=0, T=Tn, (5) 
出 口 
(a) SCP-0.743,，& =0.743 (b) SCP-0.836,， e=0.836 TT 二 L, P=0 (6) 
图 1 基本 计算 单元 i 
Fig. 1 Basic computing elements 对 称 边界 : 
au_ on_ au au 
Oy Oy Oy Oy 
Sp 
底面 加 热 面 加 以 恒定 热流 : 
oT 
Sr Ns- q (8) 
() 其 他 壁面 均 为 绝热 边界 ， 
oT 
Ss, U= uin, v=0, w=0, a7 (9) 
图 2 微 针 肋 热 沉 正面 示意 图 
Fig.2 Di f th icro-pin fin heat sink HH \、T、 
18g lagram oOo: e micro-pin Nn Neat sin 2 模拟 结果 与 讨论 
表 1 针 肋 间距 尺寸 2.1 传 热 特性 
Table 1 Dimension of pin fin spacing Nu 数 是 表示 对 流 换 热 强 度 的 无 量 纲 数 , 可 由 式 
模型 ST/mm SL/mm SD /mm (10) 得 出 : 
SCP-0.743 0.44 0.28 0.36 hd 
SCP-0.836 0.44 0.44 0.49 J (10) 
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其 中 , h 为 对 流 换 热 系数 ，W/(n2K)，d 为 当量 直 
径 , m。 

由 于 在 实际 应 用 中 ， 微 型 器 件 和 设备 的 工作 环 
境 、 条 件 和 状态 会 发 生变 化 ， 使 得 散热 器 的 热 负 荷 
发 生变 化 ， 因 此 本 文 研 究 了 不 同 热流 密度 下 微 针 肋 
热 沉 Nu 数 的 变化 情况 . 图 3 所 示 为 SCP-0.743 在 
底面 热流 密度 分 别 为 =4x105 W/m?, 8x105 W/m? 
1.2x106 W/m? 时 ，Nu 数 与 Re 数 变 化 关系 的 曲 
线 图 。 由 图 3 中 可 以 看 出 ， 相 同 热 流 密度 下 ， 随 
着 Re 数 的 增 大 ，Nwu 数 也 逐渐 增 大 ; 相同 Re 数 
下 ， 热 流 密度 的 升 高 ， 对 Nu 数 影响 不 大 . 图 4 所 
示 为 SCP-0.836 在 底面 热流 密度 分 别 为 9=4x105 
W/m?, g=8x105 WA/mm2，1.2x105 W/m? 时 ,Nu 
数 与 Re 数 变 化 关系 的 曲线 图 。 从 图 4 中 可 以 看 
出 , 在 相同 Re 数 下 ， 随 着 热流 密度 升 高 ，Nu 数 有 
较为 明显 的 上 升 。 由 此 可 见 ， 在 相同 Re 数 下 ， 对 


—s—g=4x10” W/m 
—o— g=8x105 W/m 
—A— g=1.2x10 W/m 


40 60 80 100 120 140 


Re 
图 3 SCP-0.743 在 不 同 热 流 密度 下 Nw 数 与 Re 数 关 系 
Fig. 3 Variation of Nu with Re in SCP-0.743 
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—s_4g=4x10” W/m 
—o—g=8x105 W/m 
A g=1.2x10 W/m? 


40 60 80 


100 120 140 
Re 
图 4 SCP-0.836 在 不 同 热 流 密度 下 Nw 数 与 Re 数 关 系 
Fig. 4 Variation of Nu with Re in SCP-0.836 


160 180 200 220 


不 同 空隙 率 的 微 针 肪 热 沉 ， 热 流 密度 对 热 沉 Nu 数 
的 影响 有 较 大 区 别 . 由 于 流体 流 型 对 热 沉 的 传 热 特 
性 有 很 大 影响 ， 因 此 ， 本 文 将 研究 不 同 孔 际 率 的 微 
针 肪 热 沉 中 流体 流 型 随 热 流 密度 的 变化 ， 来 分 析 流 
型 可 能 对 换 热 产生 的 影响 ， 
2.2 流动 特性 

为 方便 观察 通道 中 流体 的 流 线 图 随 Re 数 和 热 
流 密度 g 升 高 而 变化 的 趋势 ， 本 文 将 通道 沿 流 动 方 
向 ,在 入 口 、 中 段 和 出 口 处 分 别 截取 了 相同 长 度 的 
观察 区 域 . 图 5 为 入 口 雷诺 数 Reis 二 60， 热 流 密 度 
gq 二 4x105W/m? 时 , 热 沉 通道 中 流体 的 流 线 图 . 可 以 
看 出 , 在 Rein=60, 时 , 针 肪 尾部 出 现 一 对 较 小 的 旋 
涡 。 由 于 SCP-0.743 中 后 排 针 肋 距 前 排 针 肋 较 SCP- 
0.836 的 近 , 因此 , 针 肪 尾部 的 尾 涡 区 范围 较 小 , 在 针 
肪 区 出 口 处 , 由 于 SCP-0.743 的 通道 最 大 流速 (1.03 
m/s) 比 SCP-0.836 的 最 大 流速 (0.94 m/s) 大 , SCP- 
0.743 最 后 一 排 针 肋 的 尾 涡 区 明显 长 于 SCP-0.836.， 


3 


N 

(a) SCP 一 0.7 

图 5 Rein = 60, gq = 二 4x105 W/m? 时 , 通道 中 流体 的 流 线 图 
Fig. 5 Flow field in different heat sinks at 

Rein = 60, g = 4 x 105 W/m? 


随 着 Re 数 逐 渐 增 大 ， 针 肋 对 流体 的 扰动 加 剧 ， 
出 现 不 同形 态 的 涡 。 如 图 6 所 示 ， 可 以 观察 到 对 称 
涡 、 单 个 涡 、S 型 涡 等 . 不 同 空 际 率 下 ,相同 观察 区 
域内 的 涡 的 形态 也 不 同 。 

如 图 6(aj)、(b), 当 Reiw=100 时 , 在 SCP-0.743 
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中 可 以 观察 到 少量 单个 涡 、S 型 涡 的 存在 . 而 在 SCP- 
0.836 中 只 有 对 称 涡 存 在 。 随 着 Re 数 的 增 大 ， 如 图 
6(c)、(d)， 当 Rein=160 时 ，SCP-0.743 中 出 现 单 个 
涡 、S 型 涡 的 范围 增加 .。 相同 Re 数 下 , SCP-0.836 也 
开始 出 现 单个 涡 、S 型 涡 , 如 图 6(e)、 (f), 当 Rein=200 
时 , 两 种 热 沉 中 可 以 更 加 明显 观察 到 不 同形 态 的 涡 . 
出 现 这 一 现象 的 原因 是 , 对 于 SCP-0.743, 由 于 孔隙 
率 较 小 , 针 肋 排 布 较 密 , 最 大 速度 较 大 , 扰动 更 加 剧 
烈 , 相对 于 SCP-0.836, 较 小 的 Re 数 下 绕 流 针 肋 产 
生 的 扰动 已 不 能 被 务 性 力 抑制 ， 使 对 称 的 涡 演 变 为 
单个 涡 、S 型 涡 。 随 着 Re 数 的 增 大 , 通道 中 最 大 速 
度 的 增加 使 流体 的 扰动 加 剧 ， 能够 更 多 的 观察 到 针 
肋 尾 迹 区 的 涡 由 对 称 涡 变 为 单个 涡 或 S 型 涡 . 


A 
(e) 
Reis=100, 7 mmL<8.96 mm, (a) SCP—0.743, (b) SCP—0.836; 
Reis=160, 3.64 mm<L<5.6 mm, (c) SCP—0.743, (d) SCP—0.836; 
Rein=200, 3.08 mmL<5.04 mm, (e) SCP—0.743, (f) SCP 一 0.836; 
图 6 不 同 Re 数 下 流 形 的 变化 情况 : 
Fig. 6 Flow field under different Re number: 


由 于 在 实际 应 用 中 散热 右 的 热 负 茶会 发 生变 化 ， 


因此 还 需 研 究 不 同 热 流 密度 下 流体 在 通道 中 流动 时 
流 型 的 变化 。 本 文 用 数值 模拟 的 方法 ,研究 了 相同 
Re 数 时 ， 三 种 不 同 热 流 密度 下 通道 内 流体 流 型 的 
变化 情况 . 图 7 为 SCP-0.743 在 Reis 一 200 时 , 不 
同 热流 密度 下 流 型 的 变化 情况 . 对 比 图 7(a)~(c) 可 
以 看 出 ， 随 着 加 热 面 热流 密度 的 升 高 ， 针 肋 尾 部 涡 
发 展 的 速率 发 生变 化 。 以 出 口 段 为 例 ， 可 以 观察 到 
随 着 热流 密度 的 升 高 ， 尽 管 流体 条 度 减 小 ， 流 体 绕 
流 针 肋 时 黏 清 力 减弱 ， 使 壁面 效应 减弱 ， 可 以 加 强 
换 热 ， 但 由 于 前 后 排 针 肋 间距 较 小 ， 对 流体 形成 挤 
压 , 使 涡 的 出 现 相 对 沾 后 ,流体 扰动 的 剧烈 程度 下 
降 ， 削 弱 了 换 热 能 力 。 两 因素 的 共同 作用 ， 使 热流 
密度 对 Nu 数 的 影响 相对 较 小 . 图 8 为 SCP-0.836 
在 Rein=200 时 ,不同 热流 密度 下 流 型 的 变化 情况 . 
同样 以 出 口 段 为 例 , 可 以 观察 到 随 着 热流 密度 上 升 ， 
流体 温度 上 升 ， 黏 度 下 降 ， 流 体 绕 流 针 肋 时 黏 清 力 
减 小 ， 使 流体 绕 流 针 肪 时 涡 生 成 的 速率 加 快 ， 意味 
着 流体 扰动 的 加 强 , 增强 了 流体 混合 ， 加 强 了 对 流 
换 热 , 因此 Nw 数 随 热 流 密 度 升 高 而 有 相对 明显 的 
升 高 。 


(NY fg WN 人 
A \ / SN _ 
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(a) g=4x105 W/m? (b) g=8x105 W/m?(c) g=1.2x10° W/m? 
图 7 SCP-0.743 热 沉 中 流 型 随 热 流 密 度 变化 情况 : 
Fig. 7 Variation of flow field with heat flux in SCP-0.743 
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图 8 SCP-0.836 热 沉 中 流 型 随 热流 密度 变化 情况 : 
Fig. 8 Variation of flow field with heat fux in SCP-0.836 


3 结 论 


本 文 借助 CFD 模拟 软件 , 用 数值 模拟 的 方法 对 
当量 直径 为 200 wm、 和 孔 际 率 为 e =0.743 和 = =0.836 
Pea etn ee aa 

通过 分 析 模 拟 结 果 得 出 了 以 下 结 

1) 不 同 孔 隙 率 下 ， 0 
能 的 影响 不 同 . 孔隙 率 较 大 时 ， 随 着 热流 密度 上 升 ， 
热 沉 的 Nu 数 有 较为 明显 的 上 升 ; 孔隙 率 较 小 时 , 热 
流 密度 对 热 沉 Nu 数 几乎 没有 影响 ， 

2) Rein=60 时 ， 针 肋 的 尾部 可 以 形成 一 对 对 称 
的 涡 。 当 Re 数 继续 上 升 ， 流 体 绕 流 产生 的 扰动 已 
不 能 被 务 性 力 抑制 ， 在 针 肋 尾部 出 现 对 称 锅 、 单 个 


涡 、S 型 涡 等 不 同形 态 的 涡 。 随 着 Re 数 继续 增 大 ， 
单个 涡 、S 型 涡 出 现 的 更 多 . 孔 际 率 小 的 微 针 肋 热 沉 
中 更 容易 出 现 不 同形 态 的 涡 。 

3) 不 同 孔 际 率 下 ， 热流 密 度 对 流 型 的 影响 有 区 
别 . 孔隙 率 为 0.743 时 , 热流 密度 升 高 会 使 针 肋 尾部 
的 涡 滞 后 形成 ,流体 扰动 减弱 ， 抵 消 了 部 分 流体 秋 
度 降低 、 壁 面 效 应 降低 带 来 的 强化 传 热 效 果 。 孔 际 
率 为 0.836 时 , 随 着 热流 密度 升 高 , 针 肋 尾部 的 涡 加 
速 形成 ， 加 强 了 通道 内 流体 的 混合 和 扰动 ， 带 来 相 
对 较为 明显 的 强化 传 热 效 果 ,Nu 数 有 所 提升 。 
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